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ABSTRACT 

Context. The formation of water on the dust grains in the interstellar medium may proceed with hydrogen peroxide (H2O2) as an 

intermediate. Recently gas-phase H2O2 has been detected in p Oph A with an abundance of ~10""' relative to H2. 

Aims. We aim to reproduce the observed abundance of H2O2 and other species detected in p Oph A quantitatively. 

Methods. We make use of a chemical network which includes gas phase reactions as well as processes on the grains; desorption from 

the grain surface through chemical reaction is also included. We run the model for a range of physical parameters. 

Results. The abundance of H2O2 can be best reproduced at ~6xl0^ yr, which is close to the dynamical age of p Oph A. The 

abundances of other species such as H2CO, CH3OH, and O2 can be reasonably reproduced also at this time. In the early time the 

gas-phase abundance of H2O2 can be much higher than the current detected value. We predict a gas phase abundance of O2H at the 

same order of magnitude as H2O2, and an abundance of the order of 10"** for gas phase water A few other species of interest are also 

discussed. 

Conclusions. We demonstrate that H2O2 can be produced on the dust grains and released into the gas phase through non-thermal 
desorption via surface exothermic reactions. The H2O2 molecule on the grain is an important intermediate in the formation of water. 
The fact that H2O2 is over-produced in the gas phase for a range of physical conditions suggests that its destruction channel in the 
current gas phase network may be incomplete. 
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1. Introduction 



ISM: abundances - ISM: clouds - ISM: molecules - molecular processes - radio lines: ISM - stars: 



Oxyg en is the most abundant "metal" element in the cos- 
mos dSavage & SembachI Il996t lAsplund et aP l2009h . In the 
cold dense interstellar clouds, gas phase chemical models pre- 
dict that oxygen mainly resides in CO and O2 molecules 
(Herbst & Leung 1989; Millar & Herbst 1990; Wakelam et al. 
l2006h . However, although CO is ubiquitously distributed in the 
interstellar medium, O2 is not. The latter is only detected very 
recently in p Oph A at a low abundance (re lative to molecu- 
lar hydrogen) of 5x10"^ (Larsson et al. 2007), and in Orion at 
an abundance of (0.3 - 7)xl0"'' (Goldsmith et al. 201 1). On the 
other hand, the observed water (gas or ice) abundance can be 
as high as 10 (Ivan Dishoeckll20()4l) . Thus it seems that water, 
instead of O2, is a main reservoir of oxygen in addition to CO. 
When only gas phase chemistry is included, the H2 O abundance 
can be of the order 10'^ at mo st (see, for example, Bergin et alJ 
[2OOO; R oberts & Herbstll2002h for typical dai'k cloud conditions. 
The fact that O2 is over-produced and H2O is under-produced 
in gas phase chemistry suggests that adsorption onto the grain 
surfaces and the reactions on the surfaces may play an important 
role. 



On the grain surface, H2O can form through successive ad- 
ditions of hydrogen atoms to an oxygen atom: 

H H- O ^ OH, (1) 
H -H OH ^ H2O, (2) 

both of which are barrierless ( Allen & Robinsonlll977l) . It can 
also form via hydrogen addition to molecular oxygen: 



H + O2 ^ O2H, 
H + O2H ^ H2O2, 
H + H2O2 ^ H2O -I- OH. 



(3) 
(4) 
(5) 



R eaction (|3]l was assumed to have an activation barrier of 1200 K 
in lTielens & HagenI ( 1982|) . However, based on experimental re- 
;alj( 



suits, ICuppen et 1 



recently concluded that it is barrier- 
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less. Other possible formation pathways of water include the re- 
action betweenJl2_™i2Il_SP'^ route with O3 as an interme- 
diate (lTielens&Hagenlll982l) . 

In the second route described above (equations (l3]-|5]l), hy- 
drogen peroxide (HOOH, also written as H2O2, which is adopted 
in this paper) appears as an intermediate product. Thus if this 
route is indeed important, a significant amount of H2O2 might 
form on the grain, and its gas phase counterpart could also be 
detectable if effective desorption mechanisms exist. 

In the current mainstream gas phase reaction networks for 
astrochemistry, H2O2 is not efficiently formed in the gas phase. 
For example, in the 2009 version of the OSU networkQ, the only 



http://www.physics.ohio-state.edu/~eric/research_files/osu_01_2009 
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two reactions leading to the formation of H2O2 are 

H2 + O2H ^ H2O2 + H, 
OH + OH ^ H2O2 + hv. 

The first one has a large activation barrier of 10^ K, rendering it 
inactive at low temperatures. H2O2 is mainly consumed by 

H2O2 + hv ^ OH + OH, 
OH + H2O2 ^ H2O + O2H, 

the first of which is dissociation by cosmic -ray induced radia- 
tion. Other destruction channels by reacting with H and O are 
ineffective due to large activation barriers. At a temperature of 
10 K and an H2 density of 10^ cm"^, the steady-state abundance 
of H2O2 can be approximated by 

X(H202) ^ lO^X^OH) ^5x IQ-^^. 

At a higher density, the abundance of H2O2 will be even less 
because OH is less abundant in this case. With the UMIST 
RATE06 network th e abundance of H2O2 is essentially zero 
(IWoodaU et al.ll2007h . 

Thus if a substantial amount of H2O2 can be detected in the 
interstellar medium, then it must have been synthesized on the 
dust grains, rather than in the gas phase; this would also provide 
information and constraints on the formation route of H2O. 

And it was indee d recently detec t ed (for the first time) in 
the p Oph A cloud by 'Ber gman et al.l (1201 Ibh . at an abundance 
of -10 which is well above what would be predicted by the 
gas phase chemistry, indicating that chemical processes on the 
grains are responsible for this detection. Why this molecule has 
not been detected in the past seems to be a puzzle and will be 
discussed later. 

In the present work we aim at modeling the gas phase abun- 
dance of H2O2 at a physical condition relevant to p Oph A, to 
demonstrate whether the grain chemistry is able to explain its 
observed abundance. The model is also required to give consis- 
tent abundances for other species detected earlier in this region. 
Also, ice and gas-phase abundance predictions for previously 
undetected species will be done. 

The remaining part of this paper is organized as follows. In 
section |2] we describe the chemical model being used in this 
work. In section [3] we present the results of our modeling. The 
conclusions are in section |4] Appendix |A] contains an expla- 
nation to a spike-like feature in the evolution curves of some 
species. The surface reaction network we use is listed in ap- 
pendix |B] and the enthalpies of the surface species which are 
needed in the chemical desorption mechanism (see section |2]i 
are hsted in appendix ICl 

2. Chemical model 

For the gas phase chemistry, we make use of a subset of 
the UMIST RATE06 networlfl dWoodall et al.ll2007h . Species 
containing Fe, Na, Mg, CI are excluded. In total 284 gas 
phase species and 3075 gas phase reactions are included. The 
cosmic -ray ioni zation rate is taken t o be the canonical value of 
1.36x10 '^ s-i (IWoodaU et al.ll2007h . 

The surface chemical network is a combination of a selection 
of the reactions in Allen & Robinson ( 1977), Tielens & Hage^ 
(|1982|) and fHasegawa et al.l (119921) . with the rates of a few reac- 
tions updated according to the recent experimental results and/or 

^ http://udfa.net 



theoretical calculations. In total 56 surface species and 151 sur- 
face reactions are included (see appendixlBt. 

The binding energies of the surface species are either taken 
from lHasegawa & Herbstl(ll993l) . or estimated from the v alue of 
a similar species in a way similar to lGarrod et al.l (|2008|) . These 
values are applicable to bare grains, i.e., grains without an ice 
mantle. A grain will be covered by ice (typically water) as ad- 
sorption and reaction proceeds, thus these values are not always 
appropriate. Ideally, they should be varied according to the real- 
time composition of the grain. The water ice mantle mainly af- 
fects the binding energies of species with a hydrogen bond, such 
as OH and H2O. The effect of this should be minor for our pur- 
pose, because most of the reactions involving a species with a 
hydrogen bond are primarily mediated by another reaction part- 
ner which does not have a hydrogen bond. 

The barriers against surface diffusion are taken to be a fixed 
fraction of the binding energies. A range of values have been 
used for this fraction in the past, from 0.3 (JHasegawa et al. 1992|) 
through 0.5 (iGarrod et al.ll2008l) to 0.77 dRuffle & Herbst 200(|. 
We us e a value of 0.77, based on the analysis of iKatz et al] 
(IT999I) . Because our model is mainly for a relatively high tem- 
perature (~20 K) in comparison with most of the previous mod- 
els (predominantly for a temperature of ~10 K), a low diffu- 
sion barrier for the surface species would lead to an unreal- 
istic ice mantle composition. The effect of changing this pa- 
rameter is discussed later. We allow H and H2 on a dust grain 
to migrate through quantum tunneling or thermal hopping, de- 
pending on which is faster; all the heavier species are only al- 
lowed to move by thermal hopping. The quantum tunneling and 
thermal hopping rates are calculated using the formulation of 
Hasegawa et iL Ci.9-92). For calculation of the quantum tunnel- 
ing rates, we use a barrier width of 1 A. The exact value of this 
width depends on the composition and structure of the surface, 
which has not been fully quantified. 

About the activation barrier of reaction (O, as there is a big 
discrepancy between the value adopted in the past and the valu e 
proposed recently based on experiments jCuppen et alj|201^ . 
we adopt an intermediate value of 600 K. However, the effect of 
varying this parameter is also tested during the modeling, and 
will be discussed later. Reaction (|5]) has a barrier of 1400 K in 
iTielens & HagenI ( Il982ft . while in .Cuppen etalJ(l201Q) this re- 
action is barrierless. We choose to adopt the latter result in this 
case, because too high a barrier for it would result in building-up 
too much of H2O2 on the grain surface. 

Surface reactions with an activation barrier are allowed to 
proceed thermally or through quantum tunneling, depending on 
whichever is f aster. The formula used t o calculate the rates is also 
the same as in lHasegawa et alJ (Il992h . and the reaction barriers 
are assu med to have a width of 1 A, although different values are 
possible (iGarrod & Paulvll201 ih . 

The reaction rate of a two-body surface reaction 

A + B -4 C + ■■■ 

is 

[^diif(A) + ^diff(B)]/V(A)/V(B)//Vs, 

if A 5t B. Here A^s is the number of reaction sites of a grain, 
fcdiff(A) and k^iffiB) are the diffusion rates of A and B, and A^(A) 
and A^(B) are the number of species A and B on a single grain, 
respectively. If A = B, then the reaction rate should be 

^diif(A)/V(A)(/V(A) - 1)//Vs. 

With a number density of reaction sites being lO'^ cm"^, a dust 
grain with radius 0. 1 fim has a A^s of about 10^. 
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As we are mainly concerned with the gas phase abun- 
dances of several species, their desorption mechanism must be 
treated carefully, especially if they are mainly produced on the 
grains. Besides the normal thermal desorption, species can also 
get evaporated episodically when a cos mic -ray hits a grain. 
This i s treated in the same manner as in iHasegawa & Herbsn 
{Top). Furthermore, the non-thermal desorption mechanism via 
exothermic surface reactions (for brevity we call it "chemi- 
cal desorption" hereafter) proposed by Garrod et al. ( 20 0^ (see 
also IWatson & Salpetedll972t iGarrod et al.i,2007; .Cazaux et all 
is also included. Here the products of the exothermic re- 
actions on the grain have a probability to be directly ejected into 
the gas phase. The rate of such a desorption mechanism depends 
on the exoergicity of the reaction, as well as on the desorption 
energy of the products. A parameter char acterizi ng the efficiency 
of this mechanism (the parameter "a" in lGarrod et al. (2007)) is 
introduced, which we take to be 0. 1 . The yield of chemical des- 
orption is directly proportional to this "a" parameter, although it 
is not well-constrained. The value we adopt here gives a good 
match to the observational results. See section [X6] for a further 
discussion. The exoergicities of these reactions are estimated 
from the ent halpies of the reactants a nd products in the same 
manner as in I Allen & RobinsonI (1 19771) (their equations (3) and 
(4)), and the en thalpies of the specie s involv ed in these reactions 
are taken from iBinnewies & Milkd (Il999h or the NIST chem- 
istry web bookQ, or some other sporadic sources (see appendix 
0. 

However, the desorption mechanisms described above alone 
are not always sufficient to provide enough gas phase abun- 
dances for some species, especially at late times. Even if a large 
amount of a species is produced on the grain and released into 
the gas phase at early times, later it would be accreted back to the 
grain surface. If at this later time its production is no longer ac- 
tive (due to the exhaustion of the precursor species), its gas phase 
abundance cannot be maintaine d. Dust sputterin g dTielens et al.l 
11994) and photo-desorption (Oberg et al.ll2007h might help to 
release them to the gas phase, both of which should not be 
of great importance in a quiet cold dark cloud. Another possi- 
ble mechanism is that cosmic-ray induced radiation can disso- 
ciate the species on the grain, and when the fragments recom- 
bine, the products can possibly be ejected into the gas phase 
directly because of the energy release of the reaction, as de- 
scribed before. We implement this mechanism in the same way 
as in Ruffle & Herbst (200 ji) and|Garrod & Pauly (2011) (see 
also ICuppen & Herbst) l2007l) . namely, the cosmic-ray induced 
photo-dissociation rates for the surface species are taken to be 
the same as in the gas phase. Several dissociation branches from 
iGarrod et al.1 (f2008) (their Table 1) are included. 

In our model the numbers of all the species on a single grain 
are solved with the hybrid moment equati on (HME) approach 
(lDu&Parisell201l h. It has been shown in iDu & Paris3 (l201lh 
that the rate equation method can be inaccurate in some cases, 
and the HME approach provides a major improvement over the 
rate equation method. As the HME approach is relatively new, 
in several cases we also benchmarked our HME results with the 
exact Monte Carlo method (similar to the one of 'Vasvu nin et al.l 
[2009) as in Du & Parise (2011), and the agreement is satisfac- 
tory. It is impractical to run all the models with the Monte Carlo 
method because the run time would be too long. 

At present the layered structure of the grain mantle is not 
incorporated into our model. Although such a structure might 
be more realistic, and it is important for retaining some of the 



http://webbook.nist.gov/chemistry/ 



ice species, however, it is also possible that the interstellar dust 
grains may have an amorphous structure, which renders the lay- 
ered structure an inaccurate description. On the other hand, par- 
ticles landing on a grain are able to penetrate into the interior by 
sever al to tens of layers, as demonstrated by experiments (see, 
e.g., Iloppolo et al. thus although a model neglecting the 

layered structure is not accurate, one which deactivates all the 
layers below the outermost surface does not reflect the whole 
reality either In fact, for the H-addition reactions, whether the 
layered structure is taken into account or not only plays a minor 
role in determining the reaction rates in the accretion limit (i.e. 
when the accretion and evaporation processes are much slower 
than the reactions (Garrod 2008)), because the species involved 
in these reactions never build up a full layer 



3. Results and discussions 

3.1. Modeling p Oph A 

We have run the chemical model for physical parameters which 
are appropriate for p Oph A, where H2O2 is first detected by 
Ber gman et alj ( 1201 Ibi) . In Fig. [1] we show the abundances of 
several species as a function of time. In this model, a temper- 
ature of 21 K and a hydrogen density of 6x10^ cm"-' has been 
assumed, which are determined for p Oph A observationally by 
Bergman et al. (201 la b). The dust temperature and gas temper- 
ature are assumed to be the same. A fixed value of 15 for the 
visual extinction Ay has been adopted. We assume a canonical 
grain size of 0. 1 fim and a reaction site density of 10'^ cm"^. The 
dust-to-gas mass ratio is set to 0.01, and the dust grain material is 
assumed to have a mass density of 2 g cm"-'. We also assume that 
the ratio between the diffusion barrier and the binding energy is 
0.77, which is at the higher end of the values that have been used 
in the past, to give a reasonable ice composition. The initial con- 
dition is atom ic except for H2. The elemental abundances are the 
same as m Garrod & Paulvl(l2011h . 

The o bserved abundance of gas phase H2O2 is ~10""^ rela- 
tive to H2 dBergman et al.l201 lb*). This value is best matched at a 
time of ~6xl0^ yr. In the early time, before about 2x10^ years, 
the gas phase H2O2 abundance can be as high as ~5xlO"^. At 
late times the H2O2 abundance decreases to a very low value, 
due to the exhaustion of O2 on the grain and a full conversion of 
H2O2 into H2O. 

H2O2 is mainly formed through reaction (|4]l on the grain, fol- 
lowed by immediate desorption of the product into the gas phase 
caused by the reaction heat. About 7% of the produced H2O2 are 
released this way. Its gas phase abundance is determined by the 
adsorption and chemical desorption processes. The dissociation 
of gas phase H2O2 by cosmic -ray-induced photons is unimpor- 
tant in_consuming it, in c omparison with adsorption. 

Iloppolo et al.l (I2OO8I) modeled the abundance of H2O2 ice 
briefly, giving a value of lO"'"* - 10"'*' relative to molecular hy- 
drogen, depending on which energy barriers of several relevant 
reactions have been used. Our major goal is to model the gas 
phase H2O2 abundance, rather than the H2O2 ice. In our model 
results, the gas phase abundance of H2O2 is much higher than 
its ice counterpart. The H2O2 ice does not have a significant 
abundance in the later stage, being well below the upper limi t 
(5.2% with respect to H2O ice) given bv ' Boudin et alJ (Il998h . 
because it is constantly transformed into H2O by reacting with 
the accreted H atoms. However, by irradiating thin water ice film 
with low energy ions, Gomis et al. (2004) found that it is possi- 
ble to obtain an H2O2 to H2O ratio in the solid phase up to a 
few percent. This direct processing of the grain mantle by cos- 
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r.=0.1 fim, T=21.0 K, nH=6.0E + 05 cm"'' 



r.=0.1 i^m, T=21.0 K, nH=6.0E + 05 cm"'' 





Fig. 1. The time evolution of the abundances of selected species. A prefix "g" means a grain surface species (right panel), while a 
name without a prefix means a gas phase species (left panel). The dotted lines are the observed abundances (relative to H2) of O2 
(blue), H2CO (cyan), CH3OH (magenta), and H2O2 (black), respectively, in the p Oph A source. 



mic rays is not included in our model. The production of H2O2 
inside water ice i n an O2 rich environment triggered by UV ra- 
diation (IShi et al.| [201 1) should also be of little importance here. 
We notice in Fig. [1] (right panel) that in the early stage (before 
-SxlO"* yr) the H2O2 ice can achieve a rather high abundance, 
~5xlO"^ relative to H nucleus or ~10% relative to H2O. During 
this early period the water formation on the grain mainly pro- 
ceeds through reaction ^ with H2O2 as an intermediate, which 
is responsible for about half of the final water ice repository on 
the grain mantle. In the later s tage reaction (|2]) takes over In the 
results of lloppolo et al ] d2008h rtheirFig.41 we do not see a sim- 
ilar feature (i.e. a high abundance of H2O2 ice in the early stage). 
In our current model the layered structure of the grain mantle is 
not taken into account. It is possible that, if such a structure is 
considered, the inner layers with a relatively high H2O2 content 
might be maintained, which would give a value of a few percent 
for the H2O2 to H2O ratio in the solid phase. 

Methanol (CH3OH) and form aldehyde (H2CO) are a lso de- 
tected in the p Oph A SMI core dBergman et alj|201 Tal) . at an 
abundance of ~2xl0"^ and ~5xlO"^, respectively. Their abun- 
dances are also reproduced very well at a time of ~6xl0^ yr in 
our model. At early times, both CH3OH and H2CO have a high 
abundance. Their abundances also have a peak in the period be- 
tween 2x10^ yr to 10^ yr In our current network, CH3OH is 
mainly formed on the grains, and mainly through the addition of 
H atom to CH2OH, while the latter is mainly produced from the 
reaction between C and OH to form HOC followed by successive 
H additions. Thus the abundance of CH3OH decreases at very 
late times due to the depletion of atomic C (which is mainly in 
CH4 ice in the late stage). The normal formation channel through 
successive hydrogenation of CO is important at around 5x10^ - 
2x10^ yr. The gas phase H2CO mainly forms in the gas phase 
in the early stage (< 10^ yr), and mainly through the reaction 
CH3 -H O — > H2CO + H. Later it is mainly formed through suc- 
cessive hydrogenation of CO on the grain surface followed by 
chemical desorption. The abundance of methanol and formalde- 
hyde ice relative to water ice can be as high as ~20% at their 
peaks at a time of ~2xl0* yr, but falls down to a very small 
value in the late times. The late time abundances are consis- 
tent with the upper limit derived for quiescent env ironment and 
low mass young stellar objects in lGibb et al.l (l2004l) . However, in 



'Po ntoppidan et al.l (l2004t) a much higher abundance of CH3OH 
ice is observed along the line-of-sight of SVS4 (a dense cluster 
of pre-main sequence stars) which is close to a class protostar; 
this is consistent with the peak abundances in our model. 

From Fig. [1] (left panel) it can be seen that the abundance of 
gaseous O2 at an intermediate time of 6x10^ year is ~6xl0"^, 
which is within one order of magnitude of the observed abun- 
dances of 5x10"** for O2 (Larsson et al. 2007). The late-time 
abundance of O2 drops to a very low value, while its observed 
abundance is best matched at a time of ~2xl0^ yr. We notice 
that during the period ~(0.6 - 2)xl0^ yr, the abundance of O2 
ice has a prominent bump, reaching a peak abundance of ~10"^ 
relative to H2. At the same time, the gas phase O2 also reaches 
an abundance of ~(1 - 5)xlO"^. These values can be compared 
with the recent detection of gas phase O2 at an abundanc e of (0.3 
- 7)xl0 in Orion by Herschel (Gol dsmith et alj|201 Ih . Warm- 
up of the dust grain at this stage may release a large amount of 
O2 molecule into the gas phase. 

As a precursor of H2O2, O2H mainly forms from the reaction 
between O and OH on the grain, which does not have a barrier 
according to Hasegawa et al. ( 1992). The ratio between the gas 
phase O2H and H2O2 is almost constant throughout the evolu- 
tion, being approximately 3. Thus in our current network the gas 
phase O2H also has a remarkable abundance, which might be 
detectable in the future. 

Except at very early times, the grain mantle is mainly com- 
posed of water ice. The abundances of CO and CO2 are com- 
parable at an intermediate time of (0.3-l)xlO* yr, being about 
40-60% of water ice. This is in rough agr eement with the ice 
composition for interm ediate-mass YSOs (Gibb et al. 2004) (see 
also.Oberg et al]|201 lb . and is also in line with the suggestion of 
I An et al.N201lh that CO2 ice is mixed with CH3OH ice (the lat- 
ter is about 10% of the former in our model). Water ice mainly 
forms from reaction (|5]) in the early time, and from reaction (|2]) 
in the late time; the gas phase formation route of water only 
plays a minor role. On the contrary, the CO ice mantle is mainly 
from accretion of CO molecules formed in the gas phase. For 
the CO2 ice, it is mainly accreted from its gas phase counter- 
part in the early time, and its abundance is further increased 
through the reaction OH + CO — > CO2 -1- H in the late stage. 
At late times (> 3x10^ yr), most of carbon resides in the form 
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of CH4 ice, the latter being about half of the water ice. However, 
as far as we know, such a high abundance of CH4 ice (see also 
[Oarrod & Paulv 2011) has not been observed in the interstellar 
medium. 

The gas phase CO is heavily depleted, with an abundance 
~10"'' relative to H nucleus or ~1% relative to its ice counter- 
part, at an intermediate time (10^ - 10^ yr). Its abundance is 
mainly determined by the balance between the adsorption and 
cosmic -ray induced evaporation processes. At late times the CO 
ice abundance drops to a very low value. This is because CO is 
continuously hydrogenated into H2CO or CH3OH, and the dis- 
sociation of CH3OH by cosmic-ray induced photons produces 
CH3, which quickly becomes CH4 by hydrogenation. Taking 
into account the layered structure of grain mantle can retain a 
fair amount of CO ice. 

The abundance of gas phase H2O in the intermediate to late 
time is of the order 10"*^. At these times, its grain surface pro- 
duction route (|2]i followed by chemical desorption and the gas 
phase production route H30^ -1- e — > H2O + H plays a similarly 
important role. H30^ itself is mainly formed from successive 
protonation of atomic oxygen at this stage. 

The hydroxyl radical (OH) has an abundance -10 *^ - 10"^ 
in the gas phase at the intermediate to late times. These values 
are comparable to the observed abundance of ~(0.5 - l)xlO"^ 
in the envelope around the high-mass star-forming regi on W3 
IRS 5 obtained with HIFI on board Herschel ( Wampfler et alj 
I2OI lb . Although the physical condition in this source is different 
from p Oph A, however, as OH is readily produced and recycled 
in the gas phase by the reactions H30^ + e" ^ OH + 2H and 
Hj + OH H20^ + H2, grain processes will not play a dom- 
inant role in determining its abunda nce, especially when the 
temperature is not too low. However, iGoicoechea et alJ (l2006l) 
detected a much higher abundance (0.5 - l)xlO"^ of OH in 
the Orion KL outflows, which seems to require other formation 
pathways (e.g. shock destruction of H2O ice). 

A relatively high abundance (~10"^ and ~10"^) of gas phase 
and grain surface ozone (O3) is also obtained for f < 10^ yr. 
However, these values should be treated with caution, because 
the chemistry of O3 is very incomplete in our current model: 
neither the OSU gas phase network nor the UMIST RATE06 
network includes it. Possible gas phase destruction pathways in- 
volving atomic O and S (according to the NIST chemistry web 
book; see footnote O may lower its gas phase abundance sig- 
nificantly, but its ice mantle abundance should not be severely 
affected. 

The abundance of atomic hydrogen in the gas phase is quite 
high at late times as seen from Fig.[Tl which seems to be at odds 
with the usual results. In many gas phase chemical models, the 
abundance of atomic hydrogen in the gas phase is determined 
by the balance between its adsorption onto the dust grains and 
the dissociation of H2 molecules by cosmic rays. The adsorption 
process is thought to have a sticking coefficient close to unity, 
and evaporation is assumed not to occur (which is appropriate 
at low temperatures). In such a framework it can be found that 
the gas phase atomic hydrogen will always have a fixed den- 
sity of the order 1 cm"^. However, in our case with a tempera- 
ture of 21 K, evaporation is very fast and cannot be neglected. 
Furthermore, the dissociation of ice mantle by cosmic -ray in- 
duced photons generates atomic hydrogen, which enhances its 
gas phase abundance significantly in the late stage. 

We note that at a time of about 5x10"* years the abundances 
of several species change very quickly, and the abundances of 
some other species show a spike-like feature. This may resem- 
ble at first sight an erroneous behavior caused by the differential 



equation solver. So we ran the model with the same parameters 
using a Monte Carlo code (also used for benchmark purpose in 
Du & Parise (2011)) which is free from such problems, and it 
turns out that these features are genuine. A semi-quantitative ex- 
planation of this feature is in appendixlAl 

3.2. Chemical age versus dynamical time scale 

As noted before, the time of best agreement bet ween our mod- 
eling r esults and the observ at ional r esults of iBergman et al.l 
(12011 ah and iBergman et"an (1201 lb ) is at -6x10^ year 
Interest ingly, this time sca le is quite close to the time scales de- 
rived in Andre et al.l(l2007h (their Table 7). For example, the evo- 
lution time scale for p Oph A as estimated to be three times the 
free-fall time is (0.5 - 2)xl0^ years, being close to the statisti- 
cally estimated age, while the collisional time scale of 5.5x10^ 
years and the cross time scale of 8x10^ years are also of the 
same order of magnitude. 

It is important to define the starting point when talking about 
age. In the chemical evolution model described above, the whole 
system starts to evolve from a simple initial state: all the ele- 
ments except for hydrogen are in atomic form, and the grains 
are bare. How relevant is such an initial condition when we talk 
about the dynamical evolution of a cloud condensation? At a 
density of ~10"' cm"-*, a temperature of ~20 K, and a visual ex- 
tinction of 2, which are typical of diffuse or translucent molec- 
ular clouds, a chemical model starting from an atomic initial 
condition (except for H2) reaches steady-state for most of the 
species in several 10^ years. Because such a time scale is much 
shorter than the dynamical time scale of a cloud, it should not 
make much difference for a time-dependent model (i.e. one in 
which temperature and density etc. vary with time) to use either 
atomic or molecular initial conditions. For a time-independent 
cloud model (such as the one we are using) with constant phys- 
ical conditions, adopting an atomic initial condition and a high 
density is equivalent to assuming that the cloud was compressed 
from the diffuse interstellar medium very quickly. Converging 
flows in the interstellar medium might play such a role, although 
this seems unlikely in p Oph A because only very small velocity 
gradient was observed ( Andre et al...2007i) . 

However, the time of best match and the predicted abun- 
dances may also be very dependent on the value adopted for 
some of the modeling parameters that are not very well con- 
strained. It is necessary to see how the results would be changed 
if these parameters are varied. 

3.3. Effects of changing the energy barrier of the surface 
reaction H + O2 -> HO2 

The activation b a rrier of reaction (O was set to 1200 K in 
Tielens & Hagen' (1982'), which was taken from the theoretical 
calculation by Melius & Blint (1979) for the gas phase case. 
ICuppen et al.l (12010 ) concluded that this reaction has a negligi- 
ble barrier. We choose to use an intermediate value for the energy 
barrier of reaction namely, 600 K for the modeling. Here we 
study how would the uncertainties in this parameter affect the 
abundances of several species of interest. 

In Fig. |2]the abundances of several species at the best-match 
time (6x10'' yr) are plotted as a function of the activation energy 
baiTier of reaction (O. The temperature and density are fixed at 
21 K and 6x10^ cm"^, and the ratio between the diffusion energy 
barrier and the binding energy is set to 0.77. 
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The abundance of gas phase O2 is not affected by chang- 
ing the energy baiTier of reaction (O because its abundance 
is mainly determined by the gas phase production process 
O + OH — > O2 + H and adsorption. However, the abundance of 
grain surface O2 increases by five orders of magnitude as the 
barrier changes from to 1200 K, because reaction (O is one of 
the main reactions for the consumption of grain O2. However, 
even with a barrier of 1200 K for reaction (O, the abundance of 
O2 on the grain at late stage is too low to be detected. Although 
the rate coefficient of reaction O) is reduced by 6 orders of mag- 
nitude when its barrier changes from to 1200K, the abundance 
of O2H on the grain does not change significantly. The reason 
is that when reaction (O becomes slower, more O2 will build 
up on the grain, compensating for the effect of increasing the 
barrier of reaction ([3]). The abundance of gas phase H2O2 does 
not change much either, as its abundance is mainly determined 
by accretion onto the grain, and production by hydrogenation of 
O2H followed by partial chemical desorption. 




400 800 1200 

Reaction barrier of H+02 (K) 



Fig. 2. Dependence of the abundances of several species at the 
best-match time (6x10^ yr) on the reaction barrier of reaction 
(O, with density fixed to 6x10^ cm"^ and temperature fixed to 
21 K. 



3.4. Effects of changing the diffusion energy barriers 

The diffusion energy barrier of a species on the grain determines 
how fast it can migrate on the grain, so it basically determines 
the pace of the grain chemistry. Usually it is set to a fixed fraction 
(here we denote it by 77) of the binding energy of each species. 
The latter determines how fast a species evaporates into the gas 
phase. However, this parameter might depend on the material 
and morphology of the dust grain, as well as the property of the 
adsorbate itself, so it is very uncertain. Values in the range of 0.3 
- 0.77 have been adopted in the literature. We mainly used 0.77 
for our modeling. Here we investigate how different values of 77 
would affect the abundances of several species of interest. 

In Fig. [3] we plot the abundances of several species at the 
time of best-match (6x10^ yr) as a function of the ratio between 
the diffusion barrier and binding energy (//). Changing this pa- 
rameter has a large effect on the abundances of some species. 
For example, the abundance of CO2 ice is reduced with a higher 
rj. This is because at the time of concern it mainly forms through 
the reaction OH + CO — > CO2 + H, which requires the migra- 



tion of two relatively heavy species. With a low diffusion energy 
and a moderate temperature (~20 K), OH and CO can thermally 
hop quite fast, leading to a high abundance of CO2, which is not 
observed. However, if a lower temperature (~10 K) is adopted, 
the problem becomes the opposite: the mobilities of OH and 
CO are so low that it is difficult for them to meet each other 
to form enough CO2 ice, and certain intricate mechanism (e.g. 
three body reaction) h as to be introduced to account for this 
(lGarrod&Paulvll2011h . 

The abundance of H2O ice increases as 77 increases. At this 
stage it is mainly formed by hydrogenation of OH. As the mo- 
bility of atomic hydrogen on the grain is not greatly affected 
by the value of the diffusion energy barrier because it is al- 
lowed to migrate through quantum tunneling, the reaction rate 
of H -H OH — > H2O is not significantly affected by however, 
a larger rj leaves more OH available for water because a lower 
amount of it is consumed in forming CO2. This is also the rea- 
son for a higher CO (and species dependent on it such as H2CO 
and CH3OH) ice abundance when 77 is larger. The abundances of 
gas phase H2O2 and O2H also increase for a larger 77, albeit only 
mildly. 





H — h H202 


/ 


» 02H 

— 02 
gCO 

^ gC02 

— CH30H 
^ — < gCHBOH 

— « H2C0 

— H20 
gH20 





0.3 0.4 0.5 0.6 0.7 0.8 

Ratio between diffusion energy barrier and binding energy 



Fig. 3. Dependence of the abundances of several species at the 
best-match time (6x10^ yr) on the ratio between the diffusion 
and desorption energy, with density fixed to 6x10^ cm"^ and 
temperature fixed to 2 1 K. 



3.5. Dependence on the temperature and density 

The physical conditions (temperature and density) of p Oph 
A are determined by non-LTE radiative transfer modeling 
(Bergman et alj |201 lah . which is usually subject to uncertain- 
ties from many aspects, such as the excitation condition, source 
geometry, beam filling factor, etc. In this section we study how 
the uncertainties in the temperature and density of the system 
would affect the abundances of several species in our model. In 
Fig.|4]and Fig.|6]we plot the abundances of several species at the 
time of best-match (6x10^ yr) as a function of temperature and 
density. 

Apparently, temperature plays a much more drastic role than 
density. This is intuitively easy to understand because temper- 
ature enters the calculation of rates exponentially for the sur- 
face reactions. The general trend is that when the temperature 
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is either too low or too high, the grain surface chemistry tends 
to be inactive or unimportant. In the former case the mobilities 
of species other than atomic hydrogen (which migrates through 
quantum tunneling in our present model) are low, while in the 
latter case the surface abundances of many species are low due 
to elevated evaporation rates. 

As can be seen in Fig. H) the abundance of CO ice starts to 
decrease at a temperature of around 20 K. This value can be 
estimated as the temperature at which the ga s phase and grain 
surfac e abundance of CO are equal (see also iTielens & HagenI 
Il982h . when only the adso rption and evaporation p rocesses are 
taken into account (see also lHollenbach et 



or ation p n 
an i2009l) : 
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where £0 is the evaporation energy barrier of a species on the 
grain surface, v is the vibrational frequency of a species on the 
grain, Rq is the dust-to-gas number ratio, r is the grain radius, 
while m is the molecular mass of the species being considered. A 
typical value of lO'^ s"' for v, a dust grain radius of 0.1 yum, and 
a dust-to-gas ratio of 2.8x10"'^ have been adopted in deriving 
the number 60. Since the logarithmic part in this equation is usu- 
ally small for a typical gas density and temperature, the evapora- 
tion temperature can be approximated simply by Teva D - Eu/60. 
For C O, a canonical value of £0 is 1210 K (lAllen & Robinson 
11977 ). which gives an evap oration temperature of 20 K, while 
for water, an Eu of 1860 K dHasegawa & Herbsllll993h on bare 
graphite grains g i ves an evaporation temperature of ~30 K. In 
iGarrod & Herbsd (12006 ^ a much higher desorption energy of 
5700 K for water is used (appropriate for water ice), which 
gives a evaporation temperature of -95 K, close to the ob- 
served evapora tion temperature o f water in envelope surround- 
ing protostars (iMaret et al.l l2002|). The evaporation temperature 
of CH4 ice is close to that of CO ice. The evaporation time scale 
at the evaporation temperature can be estimated to be roughly 
v^' exp(£'D/revap) - lO^'^e''" s ^ 3.6 X 10^ yr. For higher tem- 
perature the evaporation will be much faster 

The abundance of CO2 ice initially increases with increas- 
ing temperature. This is because it mainly forms from the reac- 
tion CO + OH CO2 + H with a barrier of 80 K, and an in- 
crease in temperature greatly enhances the mobility of the react- 
ing species, as well as the probability of overcoming the reaction 
barrier. But when temperature increases more, the abundance of 
CO ice becomes so low that the CO2 abundance also drops (the 
evaporation temperature of CO2 is about 40 K thus evaporation 
is not responsible for the decline in CO2 ice abundance seen in 
Fig. |4]i. A similar trend is seen in other species, such as H2CO 
ice, CH3OH ice, as well as gas phase H2O2, CH3OH, H2CO, and 
H2O. On the other hand, for species efficiently produced in the 
gas phase, e.g. O2, its abundance increases with temperature due 
to the faster evaporation at a higher temperature. 

The abundances of gas phase H2O2, CH3OH, CO, and O2 
have quite a sensitive dependence on temperature at around 21 
- 24 K as seen from Fig. |4] For example, changing the temper- 
ature from 20 K to 22 K increases the abundance of H2O2 at 
a given time (6x10^ yr) by about one order of magnitude. This 
rather small change of 2 K in temperature is normally smaller 
than the accuracy of the kinetic temperature as determined from 
observational data. The dependence of the evolution curves of 
H2O2 and CH3OH on temperature can be seen more clearly in 



Fig. |5] Changing the temperature not only shifts the evolution 
curves horizontally, it also changes their shapes significantly. We 
can see that although for CH3OH it is possible to match the ob- 
served abundance at multiple stages, for H2O2 the best match is 
only possible at 3x10^ - 10^ yr (if we ignore the match in the 
very early stage). 

Regarding the density dependence. Fig. |6] shows that the 
abundances of the gas phase species at 6x10^ yr generally 
decrease with increasing density, because the accretion of 
molecules onto the dust grains is faster for a higher density. This 
does not necessarily mean that the abundances of the surface 
species always increase with density. For example, the abun- 
dances of CH4 ice and CH3OH ice decrease with higher den- 
sity, while the abundance of H2CO ice has the opposite trend. 
One important factor is the time when we look at the system. 
In Fig. Q we plot the abundances of H2CO ice and CH3OH ice 
as a function of time for different densities, while the temper- 
ature is fixed at 21 K. It can be seen that at a given time, the 
abundances of H2CO ice or CH3OH ice can either increase or 
decrease when the density is increased. We note that the evolu- 
tion curves of these species have a quasi-oscillatory feature. For 
the same species, the evolution curves have a similar shape for 
different densities, except that with a lower density the evolution 
is slower and thus the curves are shifted rightward. 

3.6. Discussions and limits of the model 

It can be seen from Fig. [1] Fig. |6] and Fig. |5] that H2O2 can 
be produced with a rather high abundance at certain evolution 
stages or with certain physical paramet ers, which is frequentl y 
higher than the current observed value (iBergman et al.ll201 Ibl) . 
especially in the early to intermediate time (~10^ - 10^ yr). One 
natural question to ask is why H2O2 has not been commonly 
detected in the interstellar medium. One possibility is that its 
spectral lines have been overlooked in the past. Another sim- 
ple explanation would be that it is over-produced in our model, 
because in the past the chemistry of H2O2 (as well as O2H 
and O3, etc.) may not have been studied in detail in the astro- 
chemistry context (especially the destruction reactions in the gas 
phase), so that its destruction routes might be incomplete. For 
example, there is no reaction in which H2O2 is destroyed by re- 
acting with H3 in the current mainstream gas phase chemical 
networks. As a rough estimate, if we take the rate of destruc- 
tion by such a reaction to be the same as the rate of the reac- 
tion H3 + H2O — > H30^ + H2, then the abundance of gas phase 
H2O2 can be reduced by about one order of magnitude. Further 
theoretical/experimental studies of the H2O2 chemistry would 
thus be very helpful, given the fact that it has been detected re- 
cently and its potential important role in the grain chemistry of 
water. A third possibility is that the rarity of H2O2 might be an 
age effect. From Fig.|5]we notice that the abundance of H2O2 is 
very high only in a relatively early stage (before ~5xlO^ yr). If 
for certain reason most of the cloud cores being observationally 
studied are older than this (due to some selection effects), then 
the H2O2 abundance in these objects would be too low to detect. 
Basically, at least three physical parameters, namely age, den- 
sity, and temperature, are relevant. A probability distribution of 
these three parameters of the cloud cores would help to give the 
detection probability of H2O2 (and any other molecules). In this 
sense, p Oph A may be considered special in the sense that it has 
a relatively high density (~10'' cm"-') and temperature (20 - 30 
K), while most dark clou ds with a hig h density (>10'* cm"^) tend 
to be very cold (<15 K) (iBergin & Taf alla 2007). An inhomoge- 
neous physical condition would make the situation more com- 
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Fig. 7. The abundances of H2CO and CH3OH ice as a function of time for different densities. The temperature is fixed to 21 K. 



plex, which may require a self-consistent dynamical-chemical 
model. However, a thorough study on these possibilities has to 
be left for future work. 

On the other hand, for CH3OH, although it has been stud- 
ied quite extensively in the past, we notice that the gas phase 
reactions associated with it contain some important (although 
not decisive) differences between the OSU09 network and 
the UMIST RATE06 network. For example, the reaction be- 
tween CH and CH3OH to form CH3 and H2CO has a rate 
2.49xl0-i°(r/300)-i ''^^ in the UMIST RATE06 network, but it 
does not exist in the OSU09 network. One possible problem is 
that the temperature of cold interstellar medium (at most several 
tens of Kelvins) is out of the indicated valid range for many re- 
actions in RATE06, and it is not clear how to extrapolate these 
reaction rates correctly, although we have closely followed the 
instructions in Woodal l et al ] (l2Q07h . In our modeling we have 
been using the RATE06 network. 

The energy barriers for the hydrogenat ion of CO and H2CO 
on the grain are both taken to be 2500 K. IWoonI (l2002h calcu- 
lated the barrier heights of these two reactions, giving a value 
of ~2740 K and 3100 K, respectively, in the case three water 
molecules are present, with zero point energy corrections added. 
If these values were adopted, then the observed abundances of 
H2CO and CH3OH can onl y be reproduced w ithin one order of 
magnitude at best. However. iGouman ^ d2011h gives a lower bar- 
rier height (~2200 K) for the hydrogenation of H2CO, which 
would give a better agreement with the observational results than 
if lWoon (2002.) were used in our model. 

The chemical desorption is very important for the abun- 
dances of the gas phase H2O2 and CH3OH. Ho wever, the ef- 
ficie ncy of this me chanism (the "a" parameter in iGarrod et al.l 
(1200 7)) is uncertain. iGarrod & Paulvl (l201ll) adopted a low value 
of 0.01 for it, to avo id over-production of some gas phase species 
dGarrod et al.ll2007h . We use a value of 0.1 in our study, and the 
abundances of H2CO and CH3OH are not over-produced, except 
possibly in the early stages of the evolution. We note that the 
temperature and density of majo r concern in our stu dy is around 
20 K and 6x10^ cm ^ while in iGarrod et all (l2007h the temper- 
ature and density are set to 10 K and 2x10^ cm"^, respectively. 
As a test we run our model with the latter physical condition, 
and in this case CH3OH and H2CO are indeed over-produced by 
about one order of magnitude. Thus it seems that the efficiency 
of chemical desorption depends on the temperature of the dust 



grain, in the sense that at higher temperature the probability that 
the product of a surface exothermic reaction gets ejected to the 
gas phase is higher However, a detailed study of this possibility 
is out of the scope of the present paper. 

Regarding the formation of water ice, iBergin et al.l (Il998h 
proposed an interesting mechanism in which water is first 
formed in the high temperature shocked gas, and then gets ad- 
sorbed onto the dust grains in the post-shock phase. This mecha- 
nism may act as an alternative or supplement to the grain chem- 
istry route. It is not our aim here to discuss to what extent this 
mechanism contributes to the water ice budget. However, we re- 
mark that even with a temperature 1000 - 2000 K, the amount of 
H2O2 produced in a pure gas phase chemistry (using the UMIST 
RATE06 network) is still much lower than the detected level. 

4. Conclusions 

With a gas-grain chemical model which properly takes into ac- 
count the desorption of grain surface species by the heat released 
by chemical reactions, we reproduce the observed abundances of 
H2O2 in p Oph A at a time of ~6xl0^ yr. The solid phase H2O2 
abundance is very low in this stage. However, a H2O2 to H2O 
ratio of a few percent might be obtained in the solid phase if the 
layered structure of grain mantle is taken into account. 

The abundances of other species such as H2CO, CH3OH, 
and O2 detected in the same object can also be reasonably re- 
produced at a time of ~6xl0^ yr. Such a time scale is consistent 
with the evolution time scale estimated through dynamical con- 
siderations. 

O2H is a precursor of H2O2 on the dust grain, and we pre- 
dict that it has a gas-phase abundance with the same order-of- 
magnitude of H2O2 and should thus be detectable. Observational 
searches for it are under way. 

For physical conditions relevant to p Oph A, water is mainly 
in solid form, being the dominant grain mantle material. Its gas 
phase abundance is only of the order of 10 according to our 
model. 

We note that the abundance of gas-phase H2O2 in our model 
results can be much higher than the current observed level for a 
range of physical conditions. This may suggest that its gas phase 
destructing channels are incomplete. Due to the potential impor- 
tant role played by H2O2 in the formation of water, its reaction 
network needs to be studied more thoroughly in the future. 
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Other uncertainties in our modeling include the ratio be- 
tween the diffusion energy barrier to the binding energy of a 
species on the grain surface, the activation energy barriers of 
certain key reactions, as well as the efficiency of the chemical 
desorption mechanism. In the present work we mainly make use 
of their canonical values, or values which give good match to the 
observational results, and we also vary them to see the effects on 
the resulting abundances, which are significant in many cases. 
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Appendix A: An explanation of the spike-like 
features in the evolution curves 

In Fig. [T] we note that at a time of ~5xlO^ yr a spike-like fea- 
ture appears in the evolution curves of some species (e.g. gas 
phase H2O2 and CH3OH), while the abundances of some other 
species (e.g. H2O2 and O2 on the grain) change very rapidly at 
the same time. While these features may appear to be caused 
by faults in the program for solving the set of differential equa- 
tions, however, after solving the same problem with the Monte 
Carlo method which is immune to such numerical instabilities, 
we find that these features are still present, indicating that they 
are genuine. How could a smooth ordinary differential equation 
system generate such an almost-singular feature? In Fig. lA.ll we 
make a zoom-in of Fig.[T](with several species added and several 
species removed). It can be seen that although the time scale of 
the spike-like feature is relatively short, the evolution is always 
smooth (except the discreteness caused by the finite sampling 
of the curve). Then what determines the appearance of such a 
feature? 

As atomic hydrogen is central to the surface chemistry, we 
first look at the most important reactions governing its abun- 
dance on the grain. The main reactions consuming atomic hy- 
drogen on the grains are (in the following a species name pre- 
ceded by a "g" means a species on the grain, otherwise it is in 
the gas phase) 

fgOH-HgOH, 

gH202, 

OH + OH, 
H2O2, 



gH + g02H 



(A.l) 



gH + gH202 



gH20 + gOH, 
H2O + OH, 



(A.2) 



while others such as the reactions with gOa and gCH20H are 
relatively unimportant in the early times. It is mainly produced 
by 

H ^ gH, (A.3) 



k4 



gOH + gCO ^ gC02 + gH, 
gH20 ^ gOH + gH. 



(A.4) 
(A.5) 
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Fig. A.l. A zoom-in plot of Fig. [T]( with several species removed 
and several added), focusing on the spike-like feature (though 
nothing looks like a spike anymore). The abundance of gH has 
been multiplied by a factor of 10'^ (to make the plot look more 
compact). 

As the abundance of atomic H in the gas phase and the abun- 
dance of gCO and gH20 change relatively smoothly, they can be 
viewed as constant in a short time scale. The main reactions for 
the consumption and production of gOH are 



gO + gOH 



k6 



fgOzH, 



gCO+gOH 



gH + gOzH 
gH + gHaOa 
gHzO- 



IO2H, 

( gC02 + gH, 
^ \C02+H, 

gOH + gOH, 

^ gH20 + gOH, 

^gH + gOH. 



(A.6) 

(A.7) 

(A.8) 
(A.9) 
(A. 10) 



From the above reaction list we may write the evolution 
equation of gH and gOH afl 

d,gU = - ^igHg02H - ^2gHgH202 

+ kiH + ;t4gOHgCO + ksgHiO, 

d.gon = - ^6gOgOH - ^7gCOgOH 

+ 2fc8gHg02H + ^sgHgHjOz + ^logHjO, 
or in a more succinct form 

d,gU =Kigll + K2gOU + bu 
5,gOH =/C3gH + /(■4gOH + b2. 



(A.ll) 



(A.12) 



(A.13) 



where 



^1 = - ^ig02H - fc2gH202, 
K2 =^4gC0, 

Ki =2fcgg02H + ^9gH202, 

K4=- k^gO - kjgCO, 
bi =/fc3H + ksgHiO, 
bi =^iogH20. 



Here for brevity we use the name of a species to denote its average 
population on a single grain; for example, if gCO=100, it mains on 
average there are 100 CO molecules on a single grain. 



If we view k\ - K4 as well as bi and b2 as constants (of course 
they are not), then equation dA.Ob can be solved exactly; the 
solution contains an exponential part and a constant part. The 
amplitude of the exponential part will be inversely proportional 
to the determinant of the coefficient matrix 



det(A:) = 



Kl K2 



As Ki - K4 are not really constant in our problem, we expect that 
when they become such that det(/c) is close to zero, a spike-like 
or jump-like behavior would appear. Namely, we require 



K1K4 

K2Ki 



(fclgOzH + ^2gH202)(^6gO + hgCO) 

(A:4gCO)(2^8g02H + fc9gH202) 

g02H/gH202 + 1 



1.2g02H/gH202 + l 



(2.5 X lO'^gO/gCO + 1) 



(A. 14) 



In the second line of the above equation the actual value of the 
parameters have been inserted. These parameters depend on the 
physical conditions. 

To satisfy this condition (at least approximately), gO/gCO 
should be very small. In Fig. lA.2l the ratio (k\K4)I{k2K2) and the 
value of 2.5xlO**gO/gCO are plotted as a function of time. The 
abundances of O, gO, g02, gH, and gOH are also plotted for 
reference (not to scale). It can be seen that the {k\K4)I{k2K2) ratio 
does decrease and approach a value of unity before the time of 
the spike/jump feature, and the gO/gCO ratio does drop to a very 
low value monotonically. 



10' 



r,=0.1 iim, T=21.0 K, nH=6.0E + 05 



cm 
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Fig. A.l. The ratio {kiK4)/(k2Kj,) and the value of 
2.5xlO'*xgO/gCO as a function of time. They are to scale. 
The abundances of several other species are plotted only for 
reference (not to scale). 



The above mathematical argument can also be understood 
intuitively. When the gO abundance is so low that reaction ( IA.6I 1 
can be neglected, gOH is only destroyed by reaction ( IA.7I ). Each 
time a gOH radical is consumed, one gH is created (if we ne- 
glect the desorption process), and this gH will quickly react with 
gH202 or g02H to create one or two gOHs. So there will be a 
net gain in the gOH abundance, leading to its fast growth, and 
the gH abundance will increase accordingly. Thus we see that 
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reaction ( IA.8b is crucial in that it produced two gOH radicals by 
consuming only one gH. 

Therefore for such a spike-like feature to occur, the abun- 
dance of gO must decrease to a low value such that the reac- 
tion between gOH and gO becomes unimportant in consuming 
gOH. Namely, we require [gO] < [gCO] ky/ke - 5 x 10"'^ The 
abundances of atomic oxygen in the gas phase and on the grain 
surface are related by O/gO - kevap(gO)/kad(0) ^ 3 x 10*", so 
equivalently we require [O] < 1.5 x 10"^ (i.e. about a factor of 
200 less than the initial O abundance) at the time of the spike- 
like feature. Atomic oxygen is mainly consumed on the grain 
surface by reacting with another gO to form gOa (for t < 10^ yr) 
or by reacting with gOH to form gOaH. Only the latter is relevant 
here. As the abundance of gOH does not change much before the 
spike-like feature, the time scale for the consumption of atomic 
oxygen can be estimated to be ln200 x (0/gO)/(A:6gOH) ^ 
10^ yr, which is of the same order of magnitude of the time of 
occurrence of the spike-like feature. 

The time scale for the endurance of the spike-like feature 
itself can be estimated to be the time scale for the exhaustion of 
g02H or gH202 (so that equations ( lA.Ul i and ( IA.12b do not hold 
anymore) by reacting with gH, which is 

l/(^igH)- lO^yr, 

where the gH population (the average number of atomic H on 
a single grain) is taken to be a median value (10"'"*) during the 
rapidly-varying period. The fact that equation (IA.14I I seems to 
hold after this period (see Fig. IA.2I ) does not mean that gH will 
keep increasing rapidly, simply because the premise of our ar- 
gument, namely equations ( lA.llI ) and ( I A. 121 ) are not a good de- 
scription of the evolution of gH and gOH anymore. 

As the abundance of H on the grain increases, almost 
all the O3 on the grain are converted into O2 and OH. O2 
molecules on the grain are then consumed by the slower reaction 
gH + g02 — > g02H, with a time scale ~10^ yr. This explains the 
prominent peak in the evolution curve of g02 (see Fig.[T]). 

Would the spike-like feature have any practical significance 
(especially observationally)? Ideally, such a short-time feature 
may be used to contrain the age of a dense cloud, by discrimi- 
ating the abundances of certain species between their early-time 
values (before the spike-like feature) and late-time values (after 
the spike-like feature). However, due to its dependence on the 
reaction network being used, which usually contains a lot of un- 
certainties and is subject to change when new experiments are 
carried out, the question whether this feature is really relevant 
for the study of interstellar medium can only be answered by 
future investigations. 
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Appendix B: The surface reaction networl< used in 
this worl( 

Refer ences. H HL92: iHasegawa et alj ( Il992h : ICetlO: lloppolo et alj 
( l2010h : CIRLIO: 'Curoen et al.' ('201(J); ICetOS: 'loppolo et alj j2008h : 
GU: [Go umans (2011); FCet09: Fuchs et al. (200^ GWH08: 
Garrod et al. (2008); AR77: Allen & Robinson (1977'); TH82: 
Tielens & Hagen (1982); GMet08: Goumans et al. (2008); ABet04: 
Atkinson et al. (2004h; RHOO: iRuffle & Herbsll ( l200(il) ; RHOl: 
Ruffle & Herbst. t20oTbl) ; 



F. Du et al.: Production of H2O2 on the grain, Online Material p 2 



Table B.l. The surface network used in this work. The photo-dissociation reactions induced by cosmic rays and the chemical desorption reactions 
are not included here. 



Num Reaction 



Branching ratio Energy barrier (K) Reference 



1 H + H^H2 1.0 

2 H + O^OH 1.0 

3 H + 02^02H 1.0 

4 H + 03^02+ OH 1.0 

5 H + 0H^H20 1.0 

6 H + O2H ^ H2O2 0.38 

7 H + O2H ^ OH + OH 0.62 

8 H + H2O2 ^ H2O + OH 1.0 

9 H + CO ^ HCO 0.5 

10 H + HCO^H2CO 1.0 

11 H + H2CO ^ CH3O 0.5 

12 H + H2CO ^ HCO + H2 0.5 

13 H + CHsO^CHaOH 1.0 

14 H + CH2OH ^ CH3OH 1.0 

15 H + HCOO ^ HCOOH 1 .0 

16 H + C^CH 1.0 

17 H + CH^CH2 1.0 

18 H-i-CH2^CH3 1.0 

19 H + CH3^CH4 1.0 

20 H + N^NH 1.0 

21 H + NH^NH2 1.0 

22 H-l-NH2^NH3 1.0 

23 H + S^HS 1.0 

24 H + HS^H2S 1.0 

25 H + H2S^HS+H2 1.0 

26 H + CS^HCS 1.0 

27 C + S^CS 1.0 

28 O + S ^ SO 1.0 

29 + 80^802 1.0 

30 + C8^0C8 1.0 

31 H + CN^HCN 1.0 

32 H + NO^HNO 1.0 

33 H + N02^HN02 1.0 

34 H + NOs^HNOs 1.0 

35 H + N2H^N2H2 1.0 

36 H + N2H2 ^ N2H + H2 1.0 

37 H + NHCO ^ NH2CO 1 .0 

38 H + NH2CO ^ NH2CHO 1.0 

39 N + HCO ^ NHCO 1.0 

40 CH + CH^C2H2 1.0 

41 + 0^02 1.0 

42 + 02^03 1.0 

43 + C0-^C02 1.0 

44 O + HCO HCOO 0.5 

45 O + HCO ^ CO2 + H 0.5 

46 O + N^NO 1.0 

47 + N0^N02 1.0 

48 + N02^N03 1.0 

49 O + CN^OCN 1.0 

50 C + N^CN 1.0 

51 N + N^N2 1.0 

52 N-i-NH^N2H 1.0 

53 H2 + OH-^H20 + H 1.0 

54 0H + C0^C02+H 1.0 

55 H + C2^C2H 1.0 

56 H + N2^N2H 1.0 

57 H + C2H^C2H2 1.0 

58 H + HOC^CHOH 1.0 

59 C + OH ^ HOC 0.5 

60 C + OH ^ CO + H 0.5 



0.0 

0.0 
600.0 
200.0 

0.0 

0.0 

0.0 

0.0 
2500.0 

0.0 
2500.0 
3000.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
860.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
650.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
1580.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
2100.0 

80.0 

0.0 
1200.0 

0.0 

0.0 

0.0 

0.0 



HHL92 

ICetlO 

Estimated 

Estimated 

ICetlO 

CIRLIO 

CIRLIO 

ICet08 

GWH08 

FCet09 

RHOO 

Gil 

FCet09 

GWH08 

AR77 

AR77 

AR77 

AR77 

AR77 

AR77 

AR77 

AR77 

HHL92 

HHL92 

TH82 

HHL92 

HHL92 

HHL92 

HHL92 

HHL92 

AR77 

AR77 

AR77 

AR77 

AR77 

HHL92 

AR77 

AR77 

AR77 

HHL92 

AR77 

ABet04 

GMet08 

GMet08 

GMet08 

AR77 

AR77 

AR77 

AR77 

AR77 

AR77 

AR77 

ABet04 

RHOl 

HHL92 

HHL92 

HHL92 

HHL92 

HHL92 

HHL92 
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Table B.l. continued. 



Num 


Reaction 


Branching ratio 


Energy barrier (K) 


Reference 


61 


CH + OH ^ CHOH 


1.0 


0.0 


HHL92 


62 


H + CHOH -> CH2OH 


1.0 


0.0 


HHL92 


63 


OH + OH H2O2 


1.0 


0.0 


HHL92 


64 


OH + CH2 -> CH2OH 


1.0 


0.0 


HHL92 


65 


C + C ^ C2 


1.0 


0.0 


HHL92 


66 


c + O2 ^ CO + 


1.0 


0.0 


HHL92 


67 


+ CH -> HCO 


1.0 


0.0 


HHL92 


68 


+ OH ^ O2H 


1.0 


0.0 


HHL92 


69 


+ CH2 H2CO 


1.0 


0.0 


HHL92 


70 


+ CH3 CH2OH 


1.0 


0.0 


HHL92 


71 


c + ^ CO 


1.0 


0.0 


HHL92 


/Z 


C + Crl — > C2rl 


i.U 


u.u 




73 


C + NH -> HNC 


1.0 


0.0 


HHL92 


74 


C + CH2 ^ C2H2 


1.0 


0.0 


HHL92 


75 


C + NH2 ^ HNC + H 


1.0 


0.0 


HHL92 


76 


N + CH ^ HCN 


1.0 


0.0 


HHL92 


77 


N + NH2 ^ N2H2 


1.0 


0.0 


HHL92 


78 


+ NH ^ HNO 


1.0 


0.0 


HHL92 
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Appendix C: The enthalpies of the surface species 

References. BM02: iBinnewies & Milkd jl999|): NIST 

Webb ook: http:// webbook.nist .go v/chem istry/; NCe tlO: iNagy et al 
(2010); K98: iKaiser et al.l ( Il998h : VS91: IVandooren et" 
1 19911) : Urll: http://chem.engr.utc.edu/webres/331f/teams- 
98/chp/Boiler%20DC/tsld007.htm 
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Table C.l. Enthalpies of the surface species considered in this work. They are used to calculate the exoergicities of the surface reactions. 



Num 


Species 


Enlhalpy (kJ/mol) 


Reference 


1 


C 


716.7 


NIST Webbook 


2 


CH 


594.1 


BM02, p238 


3 


CH2 


386.4 


BM02, p240 


4 


CH3 


145.7 


BM02, p241 


5 


CH3O 


17.0 


NIST Webbook 


6 


CH2OH 


-9.0 


NIST Webbook 


7 


CH3OH 


-201.2 


BM02, p241 


8 


CH4 


-74.9 


BM02, p241 


9 


CN 


435.1 


BM02, p247 


10 


CO 


-110.5 


BM02, p251 


11 


C02 


-393.5 


BM02, p251 


12 


cs 


280.3 


BM02, p253 


13 


H 


218.0 


BM02, p558 


14 


H2 


0.0 


BM02, p568 


15 


H2CO 


-115.9 


BM02, p240 


16 


H2O 


-241.8 


BM02, p571 


17 


H2O2 


-135.8 


BM02, p572 


18 


H2S 


-20.5 


BM02, p574 


19 


HCN 


135.1 


BM02, p239 


20 


HNC 


135.1 


Estimated 


21 


HCO 


43.5 


BM02, p239 


22 


HCOO 


-386.8 


Urll 


23 


HCOOH 


-378.6 


BM02, p240 


24 


HCS 


296.2 


K98 


25 


HNO 


99.6 


BM02, p563 


26 


HNO2 


-78.8 


BM02, p563 


27 


HNO3 


-134.3 


BM02, p563 


28 


HS 


139.3 


BM02, p567 


29 


N 


472.7 


BM02, p693 


30 


N2 


0.0 


BM02, p699 


31 


N2H 


245.2 


VS91 


32 


N2H2 


213.0 


BM02, p570 


33 


NH 


376.6 


BM02, p563 


34 


NH2 


190.4 


BM02, p570 


35 


NH2CHO 


-186.0 


NIST Webbook 


36 


NH2CO 


-13.1 


NCetlO 


37 


NH3 


-45.9 


BM02, p577 


38 


NHCO 


-101.7 


BM02, p239 


39 


NO 


90.3 


BM02, p695 


40 


NO2 


33.1 


BM02, p695 


41 


NO3 


71.1 


BM02, p695 


42 





249.2 


BM02, p733 


43 


O2 


0.0 


BM02, p741 


44 


O2H 


2.1 


BM02, p567 


45 


O3 


142.7 


BM02, p752 


46 


OCN 


159.4 


BM02, p248 


47 


OCS 


-138.4 


BM02, p251 


48 


OH 


39.0 


BM02, p566 


49 


s 


277.0 


BM02, p811 


50 


so 


5.0 


BM02, p736 


51 


S02 


-296.8 


BM02, p743 


52 


C2 


837.74 


NIST Webbook 


53 


C2H 


476.98 


NIST Webbook 


54 


C2H2 


226.73 


NIST Webbook 



